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I , Abstract 
CO 

We present a calculation of the ground state energy of massive spinor 
fields and massive scalar fields in the background of an inhomogeneous 
magnetic string with potential given by a delta function. The zeta func- 
tional regularization is used and the lowest heat kernel coefficients are 
O ■ calculated. The rest of the analytical calculation adopts the Jost function 

formalism. In the numerical part of the work the renormalized vacuum 
energy as a function of the radius R of the string is calculated and plotted 
for various values of the strength of the potential. The sign of the energy 
^ . is found to change with the radius. For both scalar and spinor fields the 

renormalized energy shows no logarithmic behaviour in the limit R — > 0, 
as was expected from the vanishing of the heat kernel coefficient A2 , which 
is not zero for other types of profiles. 



o 
o 



1 Introduction 

Recent developments in the technique for the calculation of the zero point en- 
ergy of massive fields [0 have opened interesting possibilities in the study of soft 
boundaries and smooth potentials with spherical and cylindrical geometry im- 
mersed in the vacuum of various quantum fields. The technique uses the Jost 
function of the scattering problem related to the background potential under ex- 
amination. The vacuum energy can be expressed as an integral containing the 
logarithm of a Jost function with imaginary argument: 

*e-mail: scandurr@itp.uni-leipzig.de 
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Eb = - C ^^J2 r dkik'-mlf'^^-lnm) , (1) 

here the subscript B represents the background potential, ^ and f k are respec- 
tively the integration over the momentum k and the summation over all other 
quantum numbers I; m e is the mass of the quantum field, s is a regularization 
parameter and /i a mass parameter. The Jost function fi(ik) is unique and eas- 
ily obtainable for many types of geometries with circular, spherical or cylindrical 
symmetry and for various types of potential profiles, which is the great advantage 
of this approach. 

Representation (1) for the vacuum energy still needs a renormalization. The 
zeta functional representation of (1) and its heat kernel expansion are known to be 
a very effective tool in this context (see for instance 0). The basic idea to define 
unambiguously a renormalized zero point energy, is to impose the condition that 
the vacuum energy is zero when the mass of the quantum field reaches infinity 
(seel). 

The heat kernel coefficients themselves are also of great interest, they deter- 
mine the asymptotic behaviour of the renormalized energy they are to be 
considered as an intrinsic local feature of the geometry and of the background 
under examination. 

This renormalization scheme and the above mentioned integral representation 
have been used to solve some basic configurations: the vacuum energy of scalar, 
spinor and electromagnetic fields in the background of spherical shells with hard 
and smooth potentials has been calculated [[[], [|J, j5|, |15|, ||. The investigation 
has turned recently to magnetic fields with cylindrical symmetry. In paper ][7| a 
complete analysis of a spinor field in the background of an homogeneous magnetic 
flux tube of finite radius was carried out. The vacuum energy was found to be 
negative and it did not show a minimum for any finite value of the radius. A 
natural question is if inhomogeneous magnetic fields can minimize the energy and 
render the string stable. The question was already raised in [17]. The present 
paper extends the investigation begun in |?J to an inhomogeneous magnetic string 
with delta function profile. The delta function, although not a fully realistic 
physical model, represents a simple example of inhomogeneity, which could give 
an insight into the problem of vacuum energy in magnetic backgrounds. This 
kind of "semi-transparent" was already analyzed in || for a sphere. It has some 
features in common with a smooth potential and some with a hard boundary. In 
paper || the heat kernel coefficients for a general semi-transparent boundary were 
calculated. The quantum mechanics of spinor fields in the presence of magnetic 
fluxes has been elaborated in early works f|,|| while more recent investigation 
in this direction has been motivated by the interest in the Aharonov-Bohm effect 



llj . Singular inhomogeneous magnetic fields were examined in [|12|] for the 



calculation of the fermion determinant and in [O] for the investigation of the 
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bound states of an electron, however the ground state energy was not calculated 
in those works. In our paper the ground state energy will be calculated for a 
scalar and for a spinor field. In the first part of this paper we will calculate the 
Jost functions for both fields, the heat kernel coefficients will be found and the 
energy will be renormalized imposing the vanishing of the vacuum fluctuations for 
fields of infinite mass. In the second part of the paper we will work numerically 
on the renormalized energy finding its asymptotic behaviour for small and for 
large values of the radius of the string. We will finally show some plots of the 
renormalized energy for various values of the strength of the potential. 



2 Scalar field in the background of a magnetic 
string 

2.1 Solution of the field equation 

We quantize a scalar field $ in the presence of a classical magnetic field whose 
form is that of a cylindrical shell with delta function profile. The section of the 
string is a circle with radius R. The magnetic field is given by 

B(r) = JL6(r- R) e z (2) 

where is the magnetic flux, r = \/ x 2 + y 2 and z is the axis along which the 
cylindrical shell extends to infinity. The quantum field obeys the Klein- Gordon 
equation for the scalar electrodynamics: 

(D 2 + m 2 e )<5>(x) = (3) 
where m e is the mass of the field, and 

D 2 = d 2 - V 2 - 2ieA°d t - 2ieAV - e 2 A° 2 + e 2 A 2 , (4) 

here e is the electron charge, A^ is the vector potential of the electromagnetic 
field and the convention g^ u =diag(l, — 1, — 1, — 1) is used as well as the gauge 
d^A^ = 0. We want to find a solution of eq.(3) in cylindrical coordinates z,r, (p, 
then we take the relevant operators in (4) in the following form 

^7 ( a sin( P* . o . cosip \ 

V — ► I — cos (po r c^, sm ipo r -\ c^, o z I (5) 

V 2 - \d r rd r + ^1 + dl . (6) 

The potential four vector A^ associated with the magnetic field (2) contains a 
theta-function: 



3 



27r r 

where = (— sin 0, cos0, 0). Therefore operator (4) becomes 

* - a?-(Va r + I^ + a j) + + ^L^l (8) 

here, and in the rest of this paper, f3 will represent the strength of the back- 
ground potential 

With the ansatz of the separation of the variables the scalar field is transformed 
into 

$(x) -> exp(ip t - ip z z + irrnp) ® m (p ,p z ,r) , (10) 

where p^ is the momentum four vector and m is the orbital momentum quantum 
number. Combining (4), (8) and (10) the new field equation reads 



vl-ml- V \- { - m l R))2 +-d r + d 2 r )^ m (Po,p z ,r) =0. (11) 



or 



where = v^°o — — p 2 z . The solutions of this equation are Bessel and Neu- 
mann functions. The kind of function and their coefficients are to be determined 
by means of physical considerations. We take here the regular solution which in 
general scattering theory |H| has the following asymptotic behaviour 



$ ~°J m (AT), (13) 



* ™ ^(Uk)H^(kr) + r m (k)H^(kr)), (14) 

where J m (kr) is a Bessel function of the first kind, H^_^{kr) and H^_Jkr) are 
Hankel functions of the first and of the second kind and the coefficients f m (k) 
and fm{k) are a J° s t function and its complex conjugate respectively. In the case 
of the delta function profile the regular solution reads 
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$(r) = J m (A;r)e(/2-r) + i(/ m (A;)iji 2 l /3 (A;r) + ^(A;)iji 1 l /3 (A;r))e(r- J R), (15) 
then, we can define a field <E> J in the region r < R inside the cylinder 

•i(*,r) = J m (fcr) (16) 

which is independent of the strength (3 of the potential, and a field $° in the 
region outside the cylinder rR 

= i(/ m (A;)ifi 2 l /3 (A;r) + ^(A;)iji 1 l /3 (A;r)) (17) 

which describes incoming and outgoing cylindrical waves. The conditions for the 
field at r = R will be discussed later. 

2.2 Ground state energy in terms of the Jost function and 
normalization condition 

A regularized vacuum energy can be defined as 

where the €(n,a) are the eigenvalues of the Hamiltonian operator associated with 
(4), a — ±1 being the index for the particle-anti-particle degree of freedom, while 
n includes all oder quantum numbers, s is the regularization parameter to be put 
to zero after the renormalization and \x is a mass parameter necessary to maintain 
the correct dimensions of the energy. The string is invariant under translations 
along the z axis, therefore the energy density per unit length of the string is 




where the A( n ) are the eigenvalues of the operator defined in (12) with k = 
a/Pq — ml . We perform the integration over p z in (19), getting 

^ - l^ Jrl:- 1/2) ^ (A '" ) - • <20) 

V V ' 7 (n,a) 

The next step is to transform the summation in (20) into a contour integral. We 
enclose temporarily the system in a large cylindrical quantization box imposing 
some conditions for the field at this boundary, for instance Dirichlet boundary 
conditions. Then we can express the eigenvalues A( n ) by means of the zeros of 
solution (14), which becomes an exact equation at infinity. The derivative of the 
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logarithm of the solution will then have poles at k — A( n ). We can express (20) 
through an integral whose contour encloses these poles which lie on the real axis 
k. The deformation of the contour on the imaginary axis and the dropping of the 
Minkowski space contribution allows to reach the final form. 

1 °° pea 

£SC = E / dk {k 2 + m 2 e ) 1 - s d k \nf m {ik) , (21) 

m=-oc 

where f m {ik) is the Jost function with imaginary argument and C s = (l + s(— 1 + 
2 ln(2yu)))/(27r) is a simple function of the regularization parameter. The renor- 
malization of (21) is carried out by direct subtraction of its divergent part 



esc esc esc (00\ 

c ren ~ ° °div • 

The isolation of the divergent part from the total energy will be performed via 
heat-kernel expansion as we will see in a moment. The subtracted part should 
be added in the classical part of the energy resulting in a renormalization of the 
classical parameters of the string (in paper || this procedure is well explained), 
however we do not treat the classical energy of the system here but only the 
vacuum contribution. For the analytical continuation s — > we split £^ n it into 
a "finite" and an "asymptotic" part. 



£ren ~ + ^asi (23) 

with 

1 00 f°° d 
E T=-ifs J2 J dk[k 2 -rn^- s — pn/ m (ifc)- In f™(ik)} (24) 



and 



£ Z=-2 C ' E J dk^-ml] 1 -*— In^ifc)-^, (25) 



171= — OO 



where is a portion of the uniform asymptotic expansion of the Jost function. 
The number of orders to be included in must be sufficient to let the function 



Sub = ]nf m {ik)-]nf%{ik) (26) 

fall as m~ A (or k~ A ) for k and m equally large, in this case the integral and 
the summation in (24) converge for s — > 0. To this purpose three orders in the 
asymptotics are enough. More orders would only give a quicker convergence. The 
splitting proposed in (23) immediately permits the analytical continuation s = 
in Sj c , furthermore it allows a very quick subtraction of the pole terms in the 
asymptotic part (24). £^ c s is a finite quantity. For the definition of £^ v we write 
the sum (18) as 
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1 ..2s poo 

gsc = i » / dt1 f-a/2 K f t ) (27) 

2r(s-l/2) J W 1 ; 

where if (t) is the heat kernel related to the Hamilton operator , which can be 
expanded for small t 

-tm 2 00 1 

K(t) = $>^> „ J^J^, 3 = 0,2,1... (28) 

(n) V ^ j 

The are the heat-kernel coefficients related with the Hamiltonian operator. 
By means of (27) and (28) we can expand the ground state energy in powers of 
the mass and get 

£SC = \- \? a T(s + j-2) 4 _ 2(S+J) 

L Z^ 32vr 2 r( s + i) m - 3 w 

in which the divergent contribution can be isolated 



esc m e ( 1 1 \ ?Tlg . 



e 

— ^ - + ln-^- 1 Ui +-^^3/2 
1 + (30) 



2 



327T 2 V s m 



e 



In this definition the poles are all contained in the three terms corresponding 
to the heat kernel coefficients A ,Ai,A 2 , however we included in S^ v two more 
terms in order to satisfy a normalization condition, namely that the renormalized 
ground state energy vanishes for a field of infinite mass 



lim Sr en = . (31) 

m e — >oo 

This condition fixes a unique value for the vacuum energy of a massive field. It 
is also necessary to eliminate the arbitrariness of the mass parameter \x. This 
normalization condition must of course be changed if one desires to investigate 
massless fields. In paper |15j the interested reader can find further details and 
comments about the procedure briefly summarized in this subsection. 



2.3 The Jost function and its asymptotics 

We focus our attention on the Jost function f m (k) related with the delta function 
potential of the magnetic string. 

We have defined a field in the outer region r > R and a field in the inner 
region r < R, we now need some matching conditions on the boundary r = R. 
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As the delta function is a continuous function, we will require that the field is 
continuous on the boundary. From this condition and from the field equation 
(12) it follows directly that the first derivative of the field must be continuous on 
the boundary: 

/ ®°(r)\ r=R = &(r)\ r=R 

\ d r $°(r)\ r=R = d r ¥{r)\ r=R [6Z) 
and with the use of solution (16) and (17) 



Jm{kR) = \{f m {k)H^_ p {kR) + fa(k)H<H {kR)) 



d r J m {kr)\ r=R = {f m {k)HW p {kr) + &{k)HW p {kr))\ 
Solving for f m (k) one finds 



(33) 



= 2(d r J m (kr)\ r=R H^(kR) - d r H^(kr)\ r=R J m (kR)) 
" ' ' ^ H { ^(kR)d r H^(kr)\ r=R - H^(kR)d r H^(kr)\ r=R ' 



which with the use of the Wronskian determinant for the Hankel functions JIB 
becomes 



J m (kR)H^_ p+1 - J m+1 (kR)H^_ p (kR) + JLj m ( k R)H^\{kR) 



(35) 

The corresponding Jost function on the imaginary axis can be written in terms 
of modified Bessel I and Bessel K functions 

fmiik) = i p kR [I m K m _ p+1 + I m+1 K m „p\ + iPpi m K m _i3 . (36) 

where the arguments (kR) of the Bessel functions are omitted for simplicity. Eq. 
(36) can be written in a more compact form 

f m {ik) = i p kR[l' m K m ^ + I m K' m _ p ] . (37) 

where the prime denotes the derivative with respect to the argument. This ex- 
pression holds for positive and for negative values of m. On the contrary the 
uniform asymptotic expansion, which we need for eq.(24) and (25), is a different 
function for positive or for negative m. To find it we write the Bessel I and K 
functions of eq.(36) in the form 

I m+a ((m + a)z'), K m+a ((m + a)z'), z = kR/(m + a), (38) 
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where a can be or 1 for the Bessel I function and — (3 or —f3 + 1 for the Bessel 
K function. Their expansions for large positive orders are well known for 
instance K m+a ((m + a)z') is expanded as 



fv lV ^-([i + *T 1/2 ) l 

jr[ (m + a y J 



iT m+a ((m + a)z>) ~ j 2{ J + a) ( J + ^ 2)1/4 |l 

(39) 

where 77 = \/i + z' 2 + \n(j^p^=^ ) and the are the Debye polynomials in the 
variable x. However we are interested in an expansion in powers of the variable 
m alone: an expansion of the form 

K m+a ((m + a)z') ~ (40) 

n 

where the X n are some coefficients depending on k, R and (3. Therefore we make 
the substitution z' — > z (^S^), with z = KR/m in the argument of the Bessel 
function and we rewrite its expansion as 



TS II \ A In e -(m+a)r) 

K m+a ((m + a)z) ~ J 2(m+a) fl+fz _ gs _^w* 

' ' ; , 2 , (41) 



1 ^v°° r iv "'- ([1+ ( 2 ^) ] " 1/2) 



with the obvious change for rf . Then, re-expanding in powers l/m n we get 

K m+a (kR) ~ J^exp{^ rrT n S K {n,oc,t)} , (42) 

' m n=-l 

where t = (1 + z 2 )" 1 / 2 and the functions Sxin, a,t) are given explicitly in the 
Appendix. The corresponding expansion of the Bessel I function in negative 
powers of m will be 

1 3 

I m+a {kR) == exp{ V m- n ^(n, a,t)} , (43) 

V2vrm ^1 

where the functions Sj{n,a,t) are given in the appendix. Inserting these ex- 
pansions in (36) one finds an asymptotic Jost function valid for positive m, we 
name it f^ + (ik). To find the asymptotics for negative m we must take the Jost 
function in its form 



f m (ik) = i l3 kR[I^ m K_ m+!3 _ l + I_ m _ l K_ m+l} \+ i p pI_ m K- m+li , (44) 
which is identical to (36) because of the property of the modified Bessel functions 
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h{x) = U(x) , K p (x) = I- P (x) , (45) 

where I is any natural number and p any real number. Then in (44) a large 
positive index always correspond to a large negative m and inserting (42) and 
(43) in (44) we find an asymptotic Jost function valid for negative m, we call it 

AT (*"*)• 

For m — 0, the asymptotic Jost function can be obtained from (36), using the 
expansions 
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/, M ~^(l-^ + ^-, 1 ^- 9 ' -...}, „ = 4^; (46) 
KAZ) ~ "W I ~8i~ 2!(8i)2 + • / ' (47) 



we call this contribution f^ s (ik). The finite and the asymptotic part of the energy 
defined in (24), (25) are also split into three contributions: one for positive m, one 
for negative m and one for m — 0. The positive and negative contributions can 
be summed up in a single term, but the contribution coming from m = must be 
calculated separately and summed just numerically at the end, in fact we have 



E T = -o^ E/ dk(k 2 + m 2 e )d k (]nf+(ik)- In f^{tk)) 

\m=l Jm e 

+ r dk (k 2 + m 2 e )d k (lnf (ik)-\nfS s (ik)) ,) (4* 



and 



c c = -,c. (J2 / ^(^ 2 +^) 1 - s ^in/r ± w 

\ m =l ^ 

+ f° dfc (A; 2 + m 2 ) 1 "^ In / as (^) - s£\ 

J m e / 



(49) 



where /± = f m {ik) + f- m (ik) and f%±(ik) = f™ + {ik) + f™~(ik). 

Taking the logarithm of f^ + {ik) and f°f m (ik) and re-expanding in powers of 
m we find the functions needed in the integrands of (48) (49) up to the desired 
order. We are interested in In f^ ± (ik) and ln/g s (ifc) up to the third order, they 
read 
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In/^W = ££^— , ( 51 ) 

n t 

where t = 1/(1 + (kR/m) 2 )^ and the nonzero coefficients are 

X 1A = (3 2 , X 2A = (3 2 /A , 

X 3i3 = P 2 /24 - /5 4 /12 , X 3i5 = -P 2 /2 + /3 4 /4 , (52) 

x 3J = pyiQ . 

As we mentioned above three orders are sufficient for the convergence of £j c . 
More orders could have been included in definition (50) and (51) to let the sum 
and the integral in (48) converge more rapidlyQ. 

2.4 The asymptotic part of the energy and the heat kernel 
coefficients 

Having found the Jost function related to the cylindrical delta potential an impor- 
tant part of the calculation is done. We proceed with the analytical simplification 
of £%°. The second term in (49) which we name £^ , can be quickly calculated 
inserting in it the result (50). We find 

The first term in (49), which we name here ^*s(m) contains the sum over m which 
can be carried out with the Abel-Plana formula 



, s f°° , , 1 . . f°° dm F(im) - F(-im) , N 
J>(m) = / dmF(m) - -F(0)+ / 1 ) 1 >- . (54) 



m=l 

In our case is 



F(m)= dk(k 2 + m 2 e ) 1 - s d k \nf^ s± (ik) . (55) 
Thus £as( m ) * s s P^t three addends: 



2 



OC 



O =-\c.l dm I dkiW+mlf^dkhxf^iik) (■-)()) 



£T S 2 = \C.F(0) , (57) 



1 We would like to stress that with the introduction of asymptotic expansions in our calcu- 
lation we do not approximate the vacuum energy. The total energy as defined in (22) remains 
an exact quantity. The uniform asymptotics of the Jost function is just a mathematical tool 
which permits the analytical continuation s — > 0. 
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The contributions £^ c sl an d ^2 can be calculated with the formulas given in 
appendix B. The result is 



"e 

as2 4tt/? 967rm,i? 3 48vrm pJ R 3 ' V ; 



The divergences are all contained in the first term of £^ c sl . The first term of £^ 2 
will cancel with E*° and it remains only one term containing a positive power of 
the mass. This term is the contribution to £^ v corresponding to the heat kernel 
coefficient A 3 / 2 (see eq.(30)) and therefore it will be subtracted as well as the 
pole term. 

The last contribution to £™ is £j£ 3 , whose calculation demands a little more 
work. Using the formula displayed in appendix B to calculate the integral over 
fc, we we find 



= {-l)\c s Y,X n , 3 [-m 2 - 2s r(2 - 8 )] A nd (m e R) , (61) 
where the functions A n< j(m e R) are given by 



r( s +j/2-i) 

T{j/2)xi 



rx 

Jo T 



dm mj -" 2sin[f (j-n)] 



, f°° dm m'-° 2sin[7r(l-s-ra/2)] 



(1—s—r, 



s+j/2-1 



(62) 



We have calculated these functions by partial integration for n < 3 and j < 7 
and we show them explicitly in the appendix B. By means of those functions and 
of the coefficients (52) we find 



p \ r J ( ( m 




l-e 2 



2-rvm I ffi ) " " e 



m? e — (m e Ry 



) -) -) y/ml - {m e RY 
J m J m J 
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(3 2 m e 



(63) 



16R 1 - e 27Tm - R ' 

where the prime in the integrands denotes derivative with respect to m. The last 
term in (63) goes to zero for m e — > oo and therefore must be included in the 
renormalized energy, however, given its exponential behaviour, it does not con- 
tribute to the heat kernel coefficients. The heat kernel coefficients, which we have 
calculated up to the coefficient A 7 / 2 (including four more orders in In f as± (ik)), 
read 

= 

/3 2 7T 3 / 2 

[(3 7 r-128) / 3 2 +(256-187r)/3 4 ]7T 1 / 2 

(64) 

384i? 3 

(27/3 2 -100 / 3 4 +80/3 6 )7r 3 / 2 

24576R 5 • 

We perform the subtraction proposed in (22) and we obtain the final result 



A 


= 


, A 1/2 


A 1 


= 47T/3 2 


, A 3 / 2 


A 2 


= 


, A 5 / 2 


A 3 


= 


, A 7 /2 



esc 
°as 



ff 2 



//;{//>, R) + 

\QirR 2 YItxR 2 

P 2 



+ 
+ 



24%R 2 

(3 2 m e 



h±{m e R) 
1 



(_P t_ 

\24nR 2 12-kR 2 

h 3 (m e R) 

(3 2 (3 4 



h2(m e R) 



967im e R 3 



+ 



48-irm e R 3 



16R 1 - e 2™ e ii 



The functions h n (x) are given by 

dm 



hi(x) 
h 2 (x) 
h 3 (x) 
h±{x) 



J x 1 = 



,2nm 



Vm 2 — x 2 




Vm 2 — 



(65) 



(66) 



they are all convergent integrals which can be easily numerically calculated. The 
finite part of the ground state energy given by (48) can be integrated by parts 
giving 



esc 



= ^ Si. 



dk k 



P 
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+ 



2tt 



dk k 



fo(ik) 



2kR 



(67) 



where the coefficient X n ~ are given by (52), t = 1/(1 + {kR/m) 2 )2 and 
f±(ik) = kR [iP(l' m K m - + I m K' m _ p ) + C p (r m K m +p + I m K' m+p )] , (68) 

/o(*) = i p kR[r Q K^ + I o KL ] ■ (69) 

Equations(65) and (67) are to be considered the main analytical result of this 
paper concerning the scalar field. Their sum gives the total renormalized vacuum 
energy. The sum will be performed in the numerical part of this paper. 



3 Spinor field in the background of a magnetic 
string 

3.1 Solution of the field equation 

The analysis of a spinor field in the background of a cylindrical magnetic field 
with an arbitrary profile has been performed in . The field equation for a spinor 
field 

*(r) = ( 9l( f\ ) (70) 



. 92(r) 

in the background of a translationally invariant potential, with delta function 
profile is 



P0-me d r - "-*><*-r) 



9i{r) 



d r - (A-) p(j + me J V g 2{r} 



. (71) 



The reader is referred to paper [0] for a derivation of this equation. Let us find 
the solution to (71) for one component of the spinorQ. The decoupled equation 
for the component #2 is 



where k — yp^ — rn 2 e . The solution in the region r < R is 
gZ(r) = J m (kr). (73) 

2 We are not interested here in the complete set of solutions to eq.(71), the solution for 52 0") 
is sumcient to find the Jost function of the scattering problem. The decoupled equation for the 
component gi(r) is however indispensable and it will be used later. 
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and in the region r > R 



9?(r) = \ 

here f^ m (k) and are the Jost function and its conjugate related to the 

scattering problem for the spinor field. Solutions (73) and (74) have the same 
form as those found in the scalar case for the inner and outer space. However at 
r = R the field has not the same form, as we discuss below, owing to the presence 
of the term ^5{R — r) in the field equation (72). 

The ground state energy of the spinor field in the background of the magnetic 
string is 

Eo = £ > (75) 

n,a,cr 

where the minus sign accounts for the change of the statistics, and the e( n ,a,er) are 
the eigenvalues of the Hamiltonian 

H = -i 7 V (d x i - ieMx)) + 7 °m e . (76) 

The degree of freedom a accounts for the two independent spin states. As in the 
scalar case we calculate the energy for a section of the string. The ground state 
energy density per unit length of the string in terms of the Jost function is given 
by 

oo /»oo 

gsp in = Cg J- / dk{k 2 -m 2 e y- s d k \nf s X n {ik) . (77) 

m =-00 -> m e 

The renormalization scheme is the same we introduced for the scalar case. The 
expansion of the ground state energy in powers of the mass and the definition 
of £jfj™ are the same as in (29), (30), apart from a factor —1 coming from the 
change of the statistics. The heat kernel coefficients will be of course not the 
same, we call them B n to distinguish them from those of the scalar problem. 
The normalization condition (31) remains unchanged. The ground state energy 
is split into the two parts 



(i) 



(74) 



and 



£ r=^H / dk[k*-ml]g\lnf^(ik) -In f%-^(ik)] (78) 
£ s a r = C s J2 J m dk[k 2 -ml] 1 '* -In f™-°^(ik) - £Z n , (79) 



where spin {ik) is the uniform asymptotic expansion of the Jost function taken 
up to the third order in m. In (78) the analytical continuation s — > has already 
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been performed, while in 8^ in it will be performed after the subtraction of the 
divergent portion. (79) is a finite quantity for s = 0. 

3.2 Matching conditions and Jost function 

The matching conditions for the field on the surface of the string are not the 
same as in the scalar case. More exactly the condition for the first derivative of 
the field at r = R is different from the one in the scalar case. The field is free 
inside the magnetic cylinder i.e. independent from /3, then from eq.(71) and (73) 
we have 



kg[{r) + drgiir) - ™<? 2 '(r) = 
and outside the cylinder 

kg?(r) + d r g°(r) - ?^-g°(r) = , (81) 
the continuity of the field on the boundary is required as in the scalar case: 

gi(r)\r=R = g?(r)\ r =R: g[{r)\ r=R = g° ' (r)\ r=R , (82) 
therefore, combining (80), (81) and (82) we find 



(80) 



(<9r#2( r ))lr=i? - (d r g$(r))\r=R 
Finally, the matching conditions at r = R read 



P if \ 



(83) 



gi(r)\ r=R = g°{r)\ r=R 
(d r gi(r))\ r = R - {d r g$(r) 



\r=R 



(84) 



Inserting in (84) solutions (73) and (74) we get the system 



j 1 



• J ra 2 



■C 

m-fi 



JL J 
kR Jm 



(85) 



for positive m, and the system 



JL 



Uk)Hf] m + r m {k)Hi 
Uk)Hf_ m + r m {k)H i 



j — i 



(1) 

/3—m 



kR J ~ r 



(86) 



for negative m. In (85) and (86) the argument (kR) of the Bessel and of the 
Hankel functions has been omitted for simplicity and the prime symbol indicates 
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derivative with respect to this argument. Much in the same way as we did in the 
scalar case we find the Jost functions on the imaginary axis 



fT n+ (ik) 



(87) 



i^kR [I m K m ^ 

13+1 + Im+lKm^p] , TO > 

fT n -{ik) = i- p kR[I rn K m _p +1 + I m+1 K m _p\, m<0. (88) 

The two Jost functions are identical apart form the sign on the exponent of the 
imaginary factor. However formula (77) for the ground state energy contains 
the logarithm of the Jost function and the derivative with respect to k, then the 
imaginary factor i ±/3 which is independent of k will not contribute to the energy. 

For the calculation of the asymptotic Jost function we found more convenient 
to use instead of m the expansion parameter v given by 



m + 1/2 rmfor m — 0,1,2... 
-to — 1/2 rmfor to = — 1, — 2, ... 



(89) 



with v — |, §, ... in both cases. Then the Jost functions (87) and (88) become 



f+(ik) = i^kR 
which can be expanded for large positive to 
f~(ik) = I'^kR 



m > 



TO < 



(90) 



(91) 



which can be expanded for large negative to. The asymptotic expansions of the 
Bessel I and K functions for v and k equally large are given by 



K, 



v+a±(3 



^ exp <| x n Si(n, a, t) 

' ^= ex Pi x n S K {n,a,t) \ , 
[ n _i J 



(92) 



(93) 



where x = the functions Si(n, a, t), Sx(n, a, t) are the same as in the scalar 
case and a takes the values ±| for the Bessel I function and ±| ±(3 for the Bessel 
K function. From these formulas the logarithm of the asymptotic Jost function 
can be easily calculated up to the third order and we define 



as— spin 



Yy - 

3 ' n v n ' 



(94) 



where t — 1/(1 + (kR/v) 2 )^ and the nonzero coefficients are 

n,i = P 2 , Y 2t2 = -f3 2 /A , Y 2A = (3 2 /A , 

^3,3 = P 2 /G - /? 4 /12 , ^3,5 = -7[3 2 /8 + /3 4 /4 , 

Y 3t7 = 5/5 2 /8 • 



(95) 
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3.3 The asymptotic and the finite part of the energy 

The asymptotic part of the energy can be written, using result (94), as 



roc q ° j-j 

Sir = C S J2 / dk[k 2 - rrtf-g £ Y hn - - £|f, 



(96) 



we calculate the sum over v with the help of the Abel Plana formula for half 
integer variables which can be found in the appendix. The case m = i.e. v = 1/2 
does not need to be treated separately. We have only the two contributions 



and 



r>spm 



C as2 



4tt V s V mi 



1 



P 2 m e 
16R 



3.7 



2-2s 



r(2 - s))E n;j {m e R) . 



(97) 



(98) 



The functions S n j (x) are given in the appendix. They correspond to the functions 
S nj (a;) found in paper for a generic smooth background potential. The only 
pole term is contained in and the term proportional to m e will be subtracted 
as well as the pole term, thus (97) cancels completely with the subtraction of S^ n . 
The heat kernel coefficients which we have calculated up to the coefficient B^, 
read 



B 


= 


, B 1/2 


Bi 


= 8n{3 2 


? B 3 / 2 


B 2 


= 


7 B 5 / 2 


B 3 


= 


, B 7 / 2 


B 4 


5/3 8 
1281i? 6 







-/3 2 tt 3 / 2 /R 

(3/3 2 +2/3 4 )vr 3 / 2 
64_R 3 

(135/3 2 -68/3 4 +16/3 6 )tt 3 / 2 
12288.R 5 



After the subtraction of £ d f v the asymptotic part of the energy reads 



cspm 
C as 



2{3 2 

7TR 2 



qi{m R ) + 



ft 2 



+ 



\12-jR~ ()tt/? 2 
the functions q n (x) are 



37ri? 2 6nR 2 ) 



q 3 (m e R) 



q2{m e R) 

g 4 (m e i?) 



(99) 



12nR 2 



(100) 



q x {x) 



dv 



l + e 



2tvu 



Vu 2 ^ 



x- 



dv 



l + e 



2nu 



X A 



18 




(101) 



It is interesting to note how in both scalar and spinor cases we found the final 
expression for the asymptotic part of the energy to depend only on even powers 
of p. This was to expect for physical reasons, infact, inverting the direction of 
the magnetic flux the ground state energy should not change. 

The finite part of the energy £^ pin can be hardly analytically simplified. We 
can only integrate expression (78) by parts to obtain a final form which is suitable 
for the numerical calculation: 

ef n = dkk 

77 i J m e 

where 

fi(ik) = kR{(l v+ iK, 
+ 



3,7 



t j 



(102) 



v+h-P 



)} 



(103) 



4 Numerical evaluations 

In this section we show some graphics of £ as , Sf and of the complete renormalized 
vacuum energy £ ren as a function of the radius of the string, for the scalar and 
for the spinor field. We calculate also the asymptotic behavior of S as and Sf for 
large and for small R. Since we want to study here only the dependence on R 
and on (3, we set m e — 1 for all the calculations of this section. 

4.1 Scalar field 

As a first step we rewrite in a form in which the dependence on the relevant 
parameters is more explicit: 

C c = [P 2 9i(m e R) + (3 4 g 2 (m e R)] , (104) 

with 

»iM = (-Ms) + Yi Hx) - \ h{x) + YA hi{x) -<L + 16(1 -e**) ) ' 

92W = + + ji-) . (105) 
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the asymptotic behaviour of the h n (x) functions for x — > is found to be 



h x {x) ~ -1/24 + O(x) , 

h 2 {x) ~ + i In a; - 0.0575 + , 

Tidy ^ 

/i 3 (a:) ~ - In x + 0.442 + £>(x) , 
2 



/&4(x) 



■lux + 1.826 ; 



(106) 



the corresponding behaviour of the functions g n (x) is 

gt(x) ~ 0.0317+ 0(x) , 
g 2 (x) ~ 0.0417+ 0(x) . 



(107) 



The logarithmic contributions have cancelled. This was to expect also from the 
vanishing of the heat kernel coefficient A 2 . £^ c s is proportional to R~ 2 for R — > 
and for an arbitrary (3. For R — > oo all the /i n (x) functions fall exponentially and 
so does £*°. 

The finite part £f is also proportional to R~ 2 in the limit R — > 0, it can be 
seen substituting — > fc/i? in the integrands of expression (67) 



CSC _ ' 

/ ~ 2ttR 2 



E 

rrt=l 



dkk 



2irR 2 



m e R 



dk k 



n t 



(3 2 ' 

fo(ik)\ k ^k/R ~ ^ 



(108) 



where the Jost functions are now independent of R. We name the first addend in 
(108) £-Ym an d second addend £|q, in the plots we will display them separately. 
For large R we found numerically £'j c ~ R~ 3 , which is in agreement with the 
heat kernel coefficient A 5 / 2 shown in (64), in fact the first non vanishing heat 
kernel coefficient after A 2 determines the behaviour of the renormalized energy 
for R — > oo. Below we show the plots of all the contributions to the renormalized 
ground state energy. The functions gl(x) and gl(x) are shown as well. Each 
contribution has been multiplied by R 2 so that all the curves take a finite value 
at R = 0. We found necessary to sum up to 20 in the parameter m an to integrate 
up to 1000 in the k variable in order to obtain reliable plots. All the calculations 
where performed with computer programming relying on a precision of 34 digits. 



4.2 Spinor field 

The asymptotic part of the energy is rewritten in the form 

1 



cspin 
°as 



YlllR 2 



[{3 2 e x {m e R) + (3%(m e R)} , 



(109) 
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with 



ei(ic) = (24q 1 (x) - 4g 2 (x) + 7q 3 (x) - q 4 (x)) 
e 2 {x) = 2 (g 2 (x) - q 3 (x)) . 



The asymptotic behaviour of the g„(x) functions for x 

gi(ar) 



1/48 + 0{x) , 
1 lnx + 0.635 + £>(x) 



2 

g 3 (x) ~ -lux + 1.135 + 0(x) 
2 



94 (a:) 



lnx + 3.906 + (9(x) 



is 



(110) 



(HI) 



and the corresponding behaviour of the e n (x) functions is 

ei(x) -2+ O(x) , 
e 2 (x) ~ -1 + O(x) ; 

therefore behaves like R~ 2 for small values of R. For i? 
e~ R . 

The finite part of the energy can be rewritten as 



(112) 



oo n falls like 



pspin 



nR 2 



R2 



rn e R 



3,7 



(113) 



which is proportional to R~ 2 for i? —» and to R~ 3 for i? — > oo in agreement with 
the heat kernel coefficient -85/2- We give the plots of the functions e\(x), e 2 (x) of 
£^ pm and of £^ n n , for small and for large values of the potential strength. 
The same remarks which we made in the scalar case about the numerical limits 
of summation and integration and about the precision employed in the computer 
graphics are valid here. 



5 Conclusions 

In this paper we have carried out a complete calculation of the vacuum energy 
of two different fields in the background of a magnetic string with delta function 
profile. The renormalized vacuum energy is given in terms of convergent integrals 
(eq.(65) and (67) for the scalar field, eq.(100) and (102) for the spinor field). A 
first remark can be made about the vanishing of the heat kernel coefficient A 2 
in both scalar and spinor cases. This coefficient, contributing to £di V , is n °t 
zero for a generic background potential. The vanishing of A 2 is also observed in 
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a dielectric spherical shell with a squared profile [[OJ. It could be argued that 
singular profiles possess less ultraviolet divergences than smooth profiles. This 
statement is also confirmed by the heat kernel coefficients calculated in |3|]. 

The dependence of the sign of the energy on the radius R of the string and on 
the potential strength (3 is non trivial. In the scalar case the energy is negative 
only for large values of the potential strength, while for /3 smaller than one, the 
energy shows a positive maximum (Fig. 5). In the spinor case we have almost 
an opposite situation (Fig. 10): in the region R < 1 the energy shows a negative 
minimum for (3 < 1, while it is positive for (3 > 1. However, when R becomes large 
the vacuum energy shows the same behaviour for the scalar and for the spinor 
field: it is negative for large (3 and positive for small j3. This strong dependence 
on the parameter f3 was not observed in paper |7[], where an homogeneous field 
inside the flux tube was investigated. In fact the most relevant result of our 
calculation is that the energy numerically shows a dependence on /3 4 for large (3. 
In the contributions proportional to /3 4 cancelled and the parts proportional to 
f3 2 dominated the renormalized energy. The proportionality S ren ~ (3 4 opens the 
possibility that the inhomogeneity of the magnetic field could render the vacuum 
energy larger than the classical energy of the string, when the boundary becomes 
sufficiently hard. The total energy of the system could, then, be dominated by 
the quantum contribution, in fact we have 

PA (3 2 B + {3*C 
E TOT ~ ^ + ' (H4) 

^class ^vacuum 

where A, B and C are numbers and a is the fine structure constant. However 
in the model studied here, the profile of the potential contains a delta function 
and the classical energy is infinite. It would be interesting to study an inhomoge- 
neous magnetic field which does not contain singularities in order to have a finite 
classical energy and possibly a vacuum energy depending on (3 A . Unfortunately 
with more realistic potentials the calculations become more difficult. A feasible 
model would be that of a cylindrical shell with finite thickness and with a profile 
given by a finite height box. In this case the Jost function would be expressed in 
terms of Bessel functions and hypergeometric functions. This problem is left for 
future investigation. 
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7 Appendix A: List of the functions used for the 
asymptotic expansions of the modified Bessel 
functions 

The functions Sj(n : a,t) and S k (n,a,t) used in (42) (43) and (92)are 

S/(-l,a,f) = rl + l ln {\^) > 
Sr(0,a,t) = ilnt-|ln(i±|) , 
5/(1, a, i) = — t(~ 3 + 12a 2 + 12at + 5t 2 ) , 

5/(2, a, i) = — [8a 3 t+12a 2 (-l + 2t 2 ) + a(-26t + 30t 3 ) + 3(l-6t 2 + 5t 4 )] , 
48 

5/(3, a, t) = — (((25 - 104a 2 + 16a 4 )t 3 )/3 + 16a(-7 + 4a 2 )t 4 
128 

-(531/5 - 224a 2 + 16a 4 )t 5 - (32a(-33 + 8a 2 )t 6 )/3 

-(-221 + 200a 2 )t 7 - 240at 8 - (1105t 9 )/9) ; (115) 



= -t--Iln(l^) 
Sk(0, a, t) = -hit + — In 



2 2 V 1 -^ 

S K (l,a,t) = ~^t(-3+ 12a 2 + 12at + 5t 2 ) , 

S K (2,a,t) = — [-8a 3 t + 12a 2 (-l + 2t 2 ) + a(-26t + 30t 3 ) + 3(1 - 6t 2 + 5t 4 )] , 
48 

S K (3,a,t) = — — ((25 - 104a 2 + 16a 4 )t 3 )/3 + 16a(-7 + 4a 2 )t 4 
128 

-(-531/5 + 224a 2 - 16a 4 )t 5 - (32a(-33 + 8a 2 )t 6 )/3 

-(221 - 200a 2 )t 7 - 240at 8 + (1105t 9 )/9 , (116) 

where t — 1/(1 + (kR/m) 2 ) 2 in the scalar case and t = 1/(1 + (kR/v) 2 )2 in the 
spinor case. 
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8 Appendix B: Calculation of the integrals 

The transformation of the sum in (96) into an integral has been done with the 
Abel-Plana formula for the summation over half integer numbers 



v 2 - , U r° , r° dv F(iv) - F(-iu) . x 

m=0 >/U 1/0 

The following formulas, taken from 0, have been used for the integration over 
m and fc in ££ v £ s a c s2 and in S s ^ (eq.(56),(57) and (97)) 



'aa(l) 

n— 3 s 



f. f , f2 2M-,o * J ' mr 2s r(2- s )r(i±p)r( s + ^ 

dm dk (k 2 - miY s d k — = -f- — \ ' , ' - 

Jo J me rn" 2 (Rm e r^T(j/2) 

(H8) 

where t = 1/(1 + (kR/m) 2 )? in the scalar case and t = 1/(1 + (kR/u) 2 )2 in the 
spinor case. For the calculation of the integrals over k in Eos (eq. (58)) and in 

KZ ( e q-( 98 )) we used 

from which one can derive expression (61) and the analogous expression for the 
spinor field. The functions A n j(x) and S nj (x) up to n = 3, j = 7 read 

4 f 00 dm 



Ai,i(x) = - — / - _ e2nm Vm 2 -x 2 



X 



A 2 a(x) 



71 



x 1 — e 2nx 
A 3,3(-) = ~l dm 

4 Z" 00 , 1 / / m \ ' 1 N ' 



A 3 ,s(aO = / dm - — Vm 2 -! 2 



x 2 /„, 3 V V 1 — e 2nm ) m 



4 f°° 

A 3J (x) = -— / dm — 



^(((r^)'^)'^^< 12 °> 



x 2 / T 1 + e 2 ™ 

COO / 1 1 \ ' 



4 / /ll 

^3 3(^) = o / ^ ( ^ n ) VV" 

x 2 y, V 1 + e 2 ™ 1/ 



4 f 00 1 // v 
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Figure 2: Scalar field. The three contributions to the renormalized vacuum 
energy multiplied by R 2 ■ (3~ 4 , for (3 = 0.4. The curve £|q shows the contribution 
to the finite part of the energy coming from the orbital momentum m — 0, while 
Ef m displays the contribution coming from the sum of all other m's. 
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Figure 4: Scalar field. The contributions m = 1,2. ..5 to the finite part of the 
energy multiplied by R 2 ■ f3~ 4 , for (5 = 2.2. 




Figure 5: Scalar field. The complete renormalized vacuum energy £^ n (R) 
multiplied by R 2 ■ (3~ 4 , for different values of strength of the potential. 
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Figure 6: Spinor field. The functions e 1 (x) and 62{x) contributing to the asymp- 
totic part of the ground state energy. 




Figure 7: Spinor field. The curves of the asymptotic and of the finite part of 
the energy multiplied by R 2 ■ (3~ 4 , for (3 = 0.4. 
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Figure 8: Spinor field. The curves of the asymptotic and of the finite part of the 
energy multiplied by R 2 ■ f3~ 4 , for (5 = 2.2. 
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Figure 10: Spinor field. The complete renormalized vacuum energy S spm (R) 
multiplied by R 2 ■ (3~ 4 , for different values of strength of the potential. 
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